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We construct the A resonance as a superposition of a bare A state and the ttN continuum. It is 
parametrized by three coupling constants for local 7riVA and nnNN couplings and the A mass. The 
latter incorporates the mass renormalization due to the irNA interaction, while the results depend 
only weakly, if at all, on its wave-function renormalization. Three more renormalization constants 
are needed for the derivative contact interaction. They allow one to generate the A resonance 
dynamically. A large number of fits test the quality of different model assumptions in the P33, P31 
and P13 p-wave 7riV scattering channels. 

PACS numbers: 24.30.-v, 24.10.Eq, ll.80.Gw, 11.55.-m 



The A resonance is present in most of baryonic physics, due to its low energy and its strong coupling to 
the pion-nucleon systems. [Q, ||, [| @, H> H> @] Therefore it is is crucial to have good and simple methods 
to handle the A resonance and its effects on other hadronic physics. Such method should be based on a 
local theory where the renormalization is handled properly, and divergences can be re-absorbed consistently 
in the interactions. Otherwise, the method is nothing more than a particular fit of the data, where the 
ttNA. interactions have no universal meaning. Especially this interaction should be handled properly, due 
to strength, which makes the application of perturbation theory dubious. In this paper we propose such a 
method with local interactions, which can be solved in closed form, without resorting to perturbation theory. 
The A resonance is a spin 3/2 state, which is hard to handle in field theory, as it is generally part of a non- 
renormalizable field theory, due to the high-order derivatives in the coupling and the propagator, which yield 
high-order divergences in the self-energy corrections. Furthermore, there are at least two other issues that 
blur our image of the A resonance. First, the propagator has some unknown off-shell dependence, reflected 
by a free parameter. Second, field transformations allow different representations for the interactions and 
the propagators which, however, lead to identical observable results. ||] 

In this short paper all the problems raised above are circumvented. The propagator and the interaction are 
lumped together, so there are no issues concerning field transformations that shift terms around. Further- 
more, dimensional analysis, combined with local field theory, implies that the interaction depends, in lowest 
order, on four parameters only: two TtN A couplings, one rnrNN couplings, and the A mass. The finite 
renormalization constants are fixed by the A mass and the wave function renormalization. 



The self-energy diagram || of the ttN loop correction to the A, which has an imaginary part because of the 
ttN decay, has only two types of local momentum dependence of the loop integral in the rest frame: 



due to the different terms in the traces over polarizations and spins. The attention focusses on the real 
part. The same diagram occurs in other spin-isospin channels, with, possibly, different masses and coupling 
constants, however, the momentum and energy dependence will be the same. All the singularities are 
treated a principal value integrals. The imaginary part, associated with the homogeneous solution of the 
wave equation will follow from the eigenvalue equation below. pp[ | The two different momentum dependencies 
are associated with two different vertices (see Fig. |l|). They serve as the two 7riVA couplings. Furthermore, 
there is a third coupling for the ttttNN contact interaction, which will yield the phase shift in the absence of 
the A resonance in other spin-isospin channels. In lowest order, the contact interaction is dominated by the 
pionic current {ir* jd^n] which acts on the side of the triangle in Figure |l[ The two different time-orderings 
will be associated with two different terms in the Hamiltonian. 
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FIG. 1: The four types of interactions of the Hamiltonian. The dot represents a derivative coupling. The triangle 
indicates the pion line on which the derivative acts. 

For each channel and partial wave, in the rest frame, we can write down a general Hamiltonian for low-energy 
irN scatttering based on the observations above: 

H = |A)m A (A|+ / de\irN(e))e(nN(e)\ 
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where e n k — Win^k and zero otherwise. The momentum dependence of the contact interaction is based on 
a single irirNN coupling through the pionic current ~ [it, d^ir] with the nucleon current ~ Nj^j^N, where 
the single derivative acts either on the incoming or the outgoing pion state, which leads to the asymmetric 
form. 

The coupling functions V and W, which correspond to local interactions, contain both two terms with the 
same functional form, however, different coupling constants (V n = — W n ): 
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The kinematical factors (uj^uj^)^ 1 make the interaction relativistic invariant, and the factor k 2 for n = 1 
is due to the local, derivative coupling, in the rest frame. The coupling constants g n have a dimension 
[energy 2-2 "] (GeV), and A has a dimension of [energy 1-2 ™] (GeV). All the numerical factors from angular 
integrations, etc. are lumped in with the coupling constants g n and A. The self-energy diagram, Eq. (fi]), is 
equivalent to (g n = 1): 

s{E)= r (4) 

which is again equivalent with the first order perturbation theory results for the Hamiltonian Eq. (^) . 
This single channel problem with separable, but local interactions can be solved exactly, apart from the 
divergencies for which the finite renormalization has to be determined. The coupling functions W and V 
have the same functional form, therefore the ansatz wave function for the energy lu yields: 

|*H) = a(u)\A) (5) 

de f — !— - z(u)S(u - e)J (y/W (e)/3 (u) + y/Wx(e)^(uj))\nN(e)) 

a(w)|A)+/3b(w)|w(7rJV) )+/?i(a;)K7rJV)i) . (6) 



The function z(lo) models the ratio of the homogeneous and the inhomogeneous part of the ttN scattering 
state. For this single channel scattering state it is the only observable, which is directly related to the phase 
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shift. When this wave function is inserted into the continuum eigenvalue equation = (u> — H)\^(u})) it 
yields the relations between the spectroscopic densities Pi and the occupation number a: 

= o \uj(irN)o) + pi\u(irN)i) + (w - m A )a\A) 

- (W (u) - zWo(fj))Po\A) - (WiM - zWi(w))#i|A} 

- a(|w(7rJV) ) + |w(ttJV)i)) 

- (^ H-zFoH)/3o|w(7r^)i>-(ViM-z-l/ 1 ( W ))/3 1 | W (7r7V) ) , (7) 
where the tilde refers to the Hilbert transform, without a factor of 7T : 
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Projecting on the states \uj(ttN)q) and \u>(irN)i) yields: (i, j = 0, 1) 



Pi = £i 



1 - [V (lu) - zV (u)]\Vi{u) - zVxiu)] 
inserting this back into Eq. (fjj) and projecting onto |A) yields a relation for z: 



w — TOA 



1 - [Vb(w) - *Vb(w)][Vi(w) - 



(8) 
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which yields a straightforward but lengthy quadratic equation for z. In the case of weak contact interactions 
A ~ 0, the solution can be approximated by: 
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W + W l 
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(11) 



which can be cast in the form of an analyticity corrected energy-dependent Breit-Wigner resonance. 
The general phase shift is given by: 

tan 5 = — . 



(12) 



Since there is only one channel, the scattering process is fully characterized by z through the phase shift. 
Only the Hilbert transforms of Wo = ^Vo and Wi = ^rVi still need to be determined. 
The backward or Z-diagram contributions are included, to make the results covariant and the expressions 
simpler. The singularity is treated as a principal value singularity: 
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The integral is divergent. Subtractions make the integrals finite: 
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which yield an analytical results for the Hilbert transform for an arbitrary integer value of n: 



~gw(l ~ ^2w+l, 
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where T m stands for the m-th order Taylor expansion in to 2 around lu 2 — 0, and 
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FIG. 2: The phase shift in the P33 scattering channel which contains the A resonance. The dashed-double-dotted 
line is a fit with a dynamically generated mass through the contact interaction with A. The solid line is a fit with go 
and gi. The dashed line is a fit with go, gi, and A. The long-dashed line is a fit with gi, which has a much better 
X-square value due to a better threshold behavior. The dot-dashed line is a fit with two gi couplings to two different 
resonances. 

It should noted that the real part of the function above arises through the real part of the logarithm and 
the real part of the squareroot for lu 2 > (tozv + m-ir) 2 and u> 2 < (m/y — tItt) 2 and through the product of the 
imaginary parts of both for (mjy — rn^) 2 < ui 2 < (m/y + m^) 2 . 

The subtractions correspond to two divergent constants that should be reabsorbed in the mass and the wave 
function of the bare A state: 



W a (co) 



(0) 



(17) 
(18) 

)" + Q(n+ 1). Arguably, W 



where the renormalization constants c are fixed through: W n (mA) — (u 2 — m 2 ^)" + OO r 
could have a finite wave function renormalization. However, as one can see from Eq. ( pT[ ) this renormalization 

factor C]° can be reabsorbed in an overall shift of the coupling constants g n — ► g n /c, for the only observable 



quantity z, where 



m A is approximated by 2mA (w — m^). The coupling functions W can clearly 



be identified with the bare A state; the constant part of W is indistinquishable from the mass ttia- For 
the contact interaction there is no such identification, therefore these renormalization constants will be 
determined from the data. 



III. WITHOUT BARE RESONANCE STATE 

In the case of only contact interactions and no resonance state the expressions simplify somewhat. It 
corresponds to the Hamiltonian, Eq. (|^), with Wq = and W\ — 0, which decouples the resonance from the 
scattering states. The ratio of the homogeneous to inhomogeneous solution tt/z(lu) that leads to the phase 
shift is given by: 



VoV± + V1V0 + \J (V&Vi + ViVb) 2 + - V V1) 



2V Vi 



(19) 



Depending on the coupling constant A, the phase shift can be positive or negative. In the physical range 
tojv + fn-ir < uj < 2.0 GeV the sign change in the phase shift occurs for A varying between 1.0 and 3.0. In 
that case the finite renormalization is set to zero, such that Vi vanishes at u> = 0. 

It is possible to mimic part of the resonance with a mass my through a non-trivial renormalization of Vq 
and V%: 



V = Vq + 



(20) 
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FIG. 3: The quality of the fit and the variation of the coupling constants for the P33 channel. The fit with go and g 1 
includes all the data points from threshold to the energy indicated. The circles are the \ 2 /# > the squares the values 
of go in lOOGeV, and the diamonds are the values of g\ in 2GeV. The lines only guide the eye. The xV# drops for 
large energies due to the poorer quality of the data at higher energies. It is clear from the fits that g\ dominates at 
threshold, however, cannot fit the data alone at higher energies. There are three parameters fitted: the A mass and 
the coupling constants go and gi, the smallest number of data points is nine. 

V, = VI + 4 A ' 1} + c[ XJ) lu 2 , (21) 

where m v — -(i^(ffiy) + c Q X ^) / (c^' 1 ^ lo 2 + d^V^my)). This is normally referred to as a dynamically 
generated resonance. 



IV. RESULTS 



In order to test the quality of the model a whole set of fits were performed, excluding certain parts of 
the interaction, and including others. The least-square fit is determined by minimizing the x-square value. 
However, since there is a definite non-linear dependence of the scattering data on the coupling constants 
some care is required. As expected already above, fitting with both the contact interaction and the bare 
resonance is indeterminate, as both yield a similar effect and a flat minimum in parameter space for x-square. 
However, as it turns out all the possible fits for the A resonance in the P33 channel do much worse that a fit 
with a quasi-local interaction n = \. This interaction has the wrong dimensionality, however, produces the 
best fit with two parameters: the coupling constant and the renormalized resonance mass. This is also the 
only way to yield the right threshold behavior, if the fit includes data points beyond the resonance. 
The various results are listed below and shown in Fig. ||. In the cases that the bare resonance is presence, 
i.e., go ^ or g\ ^ 0, its mass is set to zero, i.e., ttia = 0, the renormalization constants Co will account for 
the mass, with a finite renormalization. The energy for which z = and the phase shift crosses 90° is for 
every fit within a few MeV of the A resonance at 1.232 GeV. The number of data points for the fits ranges 
between 30 and 50 depending on the range. 

• For the range 1.1 GeV to 1.6 GeV, fitting only with nonzero go yields x 2 /# = 120 and go = 0.082 and 
4 0) = 1.294. 

• For the range 1.1 GeV to 1.6 GeV, fitting only with nonzero g\ yields x 2 /# = 208 and g\ = 76000. 
and e = 10100.. The results do not depend significantly on c{ \ which is subsequently set to zero. 
However, if we set — 1, we find more appropriate numbers for g\ = 3.78 and Cq 1 ' = 0.493, since for 
this particular energy d u W[ (toa) gives a finite wave- function renormalization near unity, which makes 
the energy dependence of z very weak. 

• For the range 1.1 GeV to 1.6 GeV, fitting only with nonzero A yields xV# — 140, A = 2.31, Cq A,1 ' ) = 
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FIG. 4: The quality of the fit with the pseudo-local coupling g i . The circles are the X 2 /# values. The squares are 
the values of the coupling constants in lOOGeV. 

2.26, c^' 1 ' = —1.47, and c X '°^ = 8.69. The fit yields a better threshold behavior but a worse high- 
energy behavior compared to the go + g\ fit. 

• For the range 1.1 GeV to 1.6 GeV, fitting with nonzero go and g\ yields x 2 /# = 82, go = 0.060, 
<?i = 0.79, = 0.69, and — 0.69. Both rcnormalization constants take about half the resonance 
mass, changing this fraction does not change the vf value. Changing the range of the data set changes 
the results substantially, as can be seen in Figure [| At threshold till 1.25 GeV the <?i coupling fits the 
data well. 

• For the range 1.1 GeV to 1.6 GeV, fitting only with nonzero gi yields xV# — 56, gi = 0.44, and 

(i) 

c 2 = 1.227. It is by far the best fit archicved. 

• For the extended range 1.1 GeV to 1.8 GeV, the fit even improves, due to poorer quality data: x 2 /# = 

(-) 

45 and gi = 0.44, and c 2 — 1.227. The actual values of the coupling constant and the renormalization 
constant do not change. For the whole range the two parameters gi and fit the data with a 
consistent coupling strength of gi — 47 ± 2, as can be seen in Figure ||. Fits with data below the 
inelastic threshold the x 2 /# ranges between 16 and 23. 

• A further dramatic improvement is archieved with an additional resonance. For the range 1.1 GeV to 
2.0 GeV, the fit yields: x7# = 13, g a i = 0.48, c Q a '^ = 1.231, g b i = 1.75, and c£'^ = 3.74. The 
second resonance is around 2.5-4.0 GeV, however, its position cannot be determined accurately from 
the data, since it lies outside the range. 

Although, consistency with non-relativistic treatment, where the power n in W n is directly related to the 
partial wave, of the A resonance would require that go is zero, it is impossible to fit the P33 data properly 
with that assumption. Attempts to make sensible fits with g\ and two bare states failed. The accurate 
data at low energies forces a gi interaction in the case of two bare states. However, the P31 and the P13 
channels do give results consistent with go — 0. A number of fits were made, which are listed below, and 
show respectively in Fig. || and Fig. [| 

• For the P31 a fit with g\ only gives excellent results from the threshold up to the first structure in the 
data at 1.55 GeV. The x-square is comparable to unity. 

• Similarly for the P31 channel, a fit with A only gives good results too. However, the second solution of 
the quadratic equation starts to interfere with the physical solution at 1.6 GeV. 

• Combining g\ and A interactions does not give a significant improvement with respect to either fit. 
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FIG. 5: The phase shift in the P31 scattering channel. The long-dashed line is a fit with g\ only. The dot-dashed line 
is a fit with A only. The solid line is a fit with go, gi, and A, which does not yield a significant improvement at low 
energies. 

• For the Pi 3 the data is rather coarse. Some structure appears at 1.4 GeV. Both g\ on its own, as A on 
its own give a fit with a x-square comparable to unity. The shape of the fit depends on the range of 
the data points, however, as the x-square is comparable to unity, no conclusions can be drawn. 

• Combining both fits for the Pi 3 data leads to an improvement as some of the structure at 1.7 GeV is 
reproduced, however, as the threshold behavior fails, it is considered an artificial improvement. 

Since in the case of the P31 and the Pi 3 channels there is no resonance that fixed a scale, there is considerable 
freedom to scale the coupling constants and renormalization constants and still produce the same results. 
However, especially for a non-zero A, there seems to be no obvious rebundant parameters in the model, wich 
can be fixed a priori. The different coupling constants and renormalization constants are not quoted as they, 
on their own, have lost their significance. 

V. CONCLUSIONS 

A local field theory without cut-offs or form factors should fit the data at all energies. Deviations can only 
occur if other decay channels open, or higher-order terms in the interaction are important. Therefore we 
expect the data to be fitted well with our model below the two-pion production threshold, which is the case 
for the P31 and the Pi 3 channels. In the same theory resonance states can be pushed out of the physical 
energy region, and resonances can be generated dynamically through a contact interaction. Eventually it is 
the quality of the fit that determines which is the likely process. For the P33 channel the bare state gives a 
better results, for the P31 and the Pi 3 channels both the contact interaction and the bare resonance state 
whose pole is pushed away from the real axis through a strong interaaction do a good job fitting the data 
below 1.5 GeV. 

From the results it is clear that one cannot simply add a resonance structure to a contact interaction 
incoherently. The effects of the contact interaction are not like what one would expect from perturbation 
theory. When the interaction is summed to all orders, the resonance dominate the results. Any other, 
perturbative, effects will be irrelevant for the phase shift around the resonance peak. In all the different 
versions of the model, where some interactions are switched off, it is impossible to create the effect of a rising 
phase shift beyond 180° after the A resonance, usually associated with a background, without distorting the 
resonance itself. It can only consistently be explained by a second resonance, which is found around 2.5-4.0 
GeV. 

Furthermore, from the data the Wi coupling function prevails over the Wo and the W\ coupling functions, 
or even a combination of both with five fitting parameters, compared to the two of the Wi coupling alone. 
The same model, with different values for the coupling constants and the renormalization constants, fit the 
low-energy Pi 3 and P31 data well, which have no resonance structure. The fits give a model in with the 
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FIG. 6: The phase shift in the P13 scattering channel. The dashed line is a fit with gi only. The long-dashed line is 
a fit with A only. The solid line is a fit with go, gi, A. 



derivative coupling dominates, associated with the g\ coupling constant. 
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